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Abstract 

A simple recursive scheme for parameterisation of n-by-n unitary matrices is presented. 
The n-by-n matrix is expressed as a product containing the (n-l)-by-(n-l) matrix and a 
<^ unitary matrix that contains the additional parameters needed to go from n-1 to n. The 

procedure is repeated to obtain recursion formulas for n-by-n unitary matrices. 

<n : 

^ ; 1 The parameterisation 

o\ : 

Tj- , It has been known since a long time that unitary transformations play a central role in physics. An 
^ \ excellent example is Wigner's paper of 1939 pQ which has had a great impact in the development 
O ' °f physics and is still important. 

It is also a known fact that a general n-by-n unitary matrix X^ 1 ' may be expressed as a product 
of three unitary matrices, 
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x (n) = (i) 



■ where the matrices $ are diagonal unitary matrices, 
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is defined analogously; the a's and /5's being real. The matrix has n 2 real parameters. 
In the following, for simplicity, the word parameter stands for real parameter. The quantities a 
and (3 take care of 2n — 1 parameters of X^ because only the sums + j3j enter, where i and j 
run from 1 to n. The remaining (n — l) 2 parameters reside in the non-trivial matrix which 
is the subject of this study. 

We start by putting = 1, whereby X^ 1 ' = e^ ai+/31 ) is the most general one-by-one unitary 
"matrix". For n > 2, we write the matrix V^ n ' in the form 

v(n) _ ( + (1 - c n )\A^ >< B^)\ s n \A^ > \ . 

V " s n <B^)\ c n {6) 
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Here we have introduced an angle denoted by 6 n and have used the common notation s n = sin8 n , 
c n = cos9 n . The complex vectors "A^™ -1 )" and have each n-1 components, i.e., 
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Furthermore, 
and 



A^-V and S^" 1 ) are not arbitrary but are required to satisfy the conditions 



< y 4( n ~ 1 )|yl( ri ~ 1 ) >= i ; \B^ n ^ >= — \/( rt_1 )t| 7 4( ri - 1 ) > 



(5) 
(6) 

(7) 



whereby 

< B^-^lB^-V >= 1, \A {n - l) >= -^(™- 1 )| J B ( "- 1) > (8) 

We can easily check that if the matrix V^ n ~ l \ in Eq.fjSJ), is unitary so is V^ n \ In order for 
to be the most general n-by-n unitary matrix, modulus the phase matrices $, it must have 
the required number of parameters. The vector A^™ -1 ), having n — 1 complex components, would 
seem to represent 2(n — 1) parameters. But it has only 2(n — 2) because it is normalised and its 
overall phase can be absorbed into the matrices $, i.e., the transformation 



yields 
and 



\ A (n-i) > ^ e ^(«-i) > 



(9) 
(10) 
fill 



V (n) -> $(0, 0, e^") V (n) $(0, 0, e in ) 

The parameter counting, therefore, goes as follows. On the LHS of Eq.© we need to have (n — l) 2 
parameters. On the RHS, we have (n — 2) 2 from V^ n ~^ and 2{n — 2) from the vector A^ n ~^ Thus, 
together with the angle 9 n , the number of parameters is (n — 2) 2 + 2(n — 2) + 1 which equals 
(n — l) 2 as required. 

We may use the relation (J7J) between A^ 1 ^ and B^ 1 ^ to rewrite the matrix in terms of 
either A^ 1 ^ or B^ 1 ^. In terms of A^ n_1 \ we have 
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1 - (1 - c„))|^ (n_1) >< A^\ s n \A^ > 
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While writing the matrix in terms of B^ n ^ yields 



y(n-i) o U 1 - (1 - c n ))\B^ n -^ >< B^- 1 ^ - s n \B( n -^ > 



1 J \ s n < B( n ~V 



1 



(13) 



These relations allow for a systematic construction of unitary matrices order by order. By repeating 
the above procedure for the matrix y( n_1 ) in terms of A^ n ~ 2 ^ and B^ n ~ 2 \ and following down the 
chain we find the recursion formulas that we are looking for, 

V {n) = K,n-xK,n-1-A n , 2 K,X (14) 

y( n ) = B njl B ni 2..-B n)n _ 2 B n>n -\ (15) 

The matrices A n ^ n _\ and B n>n -\ were previously defined in Eqs. and (fT3|l . For j < n — 1 we 
have 

l-(l-c 3+1 ))|^)><#| s i+ i|^)>\ 
A nJ = ( V -« i+1 <i4W| c J+1 J | (16) 

I n -j-\ 

\-{l-c j+x ))\B^><B^\ -s m \B®> \ 
B nd = ( V *i+i<£W| J ) (17) 

Here I n -j-i is the unit matrix of order n — j — 1. The two unitary matrices A n .j and £? nj - are 
related by 

f o \ / y&')t 
2 Simple examples 

The simplest case is n = 2 for which we take \A^ >= 1 whereby \B^' >= — 1. Using y^ = 1 
we obtain, from Eqs. fT2|) and fT3|) 

(19) 



Vt2) = ^ = B » = ( -I I ) 



This is the familiar rotation matrix ^2(^2), #2 being the rotation angle in two dimensions. 

The next simplest case is n = 3 for which we may either use the "mixed form" or the pure 
forms. For the mixed form we have 

y(3) = ( R2(0) + (1 - C 3 )\AW >< S 3 \A {2) > \ (20) 



S3<5 (2) | 

Here R2{9) is again the rotation matrix in Eq. (jl9j) and we may put 



IAV >= ( M , |B« >= ( I ) (2!) 



From Eq.fJZJ) follows that 

|A (2) >= -R 2 (6)\B {2) >, \B® >= -R 2 (-e)\A® > (22) 
Hence these vectors represent two parameters, for example 

^ > y sin^e lS j 
where 7 and 5 are real. The pure forms are also obtained very simply, for example 

^Sr"' 21 ' S3| C > ) ( R t ] \) <*> 

where is as defined in Eq. ()23|) . Evidently, depending on the application one has in mind some 
choices may be more convenient than others. This is demonstrated in Refs.j2] and jSj which deal 
with the so called quark and lepton mixing matrices. The essential point is that is described 
in a rather simple fashion by four parameters as it should be. For the case of n = 4 we could, for 
example, take A^ 3 ' to be 

(cosp \ 
sinpcosae iSl (25) 
sinpsinae %&2 J 

where p, a, Si and 5 2 are the four parameters needed to define the most general A^ 3 \ 

In principle, the above recursive procedure may "easily" be extended to much larger n with 
the help of computers. 
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